Spin dephasing in n-typed GaAs quantum wells 
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We perform a many-body study of the spin dephasing due to the D'yakonov-Perel' effect in n- 
typed GaAs (100) quantum wells for high temperatures (> 120 K) under moderate magnetic fields 
in the Voigt configuration by constructing and numerically solving the kinetic Bloch equations. We 
include all the spin conserving scattering such as the electron-phonon, the electron-nonmagnetic 
impurity as well as the electron-electron Coulomb scattering in our theory and investigate how the 
spin dephasing rate is affected by the initial spin polarization, temperature, impurity, magnetic field 
as well as the electron density. The dephasing obtained from our theory contains not only that due 
to the effective spin-flip scattering first proposed by D'yakonov anc | Perel' [Zh. Eksp. Teor. Fiz. 60, 
1954(1971) [Sov. Phys.-JETP 38, 1053(1971)]], but also the recently proposed many-body dephasing 
due to the inhomogeneous broadening provided by the DP term [Wu, J. Supercond.dncorp. Novel 
Mechanism 14, 245 (2001); Wu and Ning, Eur. Phys. J. B 18, 373 (2000)]. We show that for 
the electron densities we study, the spin dephasing rate is dominated by the many-body effect. 
Equally remarkable is that we are now able to investigate the spin dephasing with extra large spin 
polarization (up to 100 %) which has not been discussed both theoretically and experimentally. We 
find a dramatic decrease of the spin dephasing rate for large spin polarizations. The spin dephasing 
time (SDT), which is defined as the inverse of the spin dephasing rate, we get at low initial spin 
polarization is in agreement with the experiment both qualitatively and quantitatively. 

PACS numbers: PACS: 71.10.-w, 67.57.Lm, 72.25.Rb, 73.61.Ey 



I. INTRODUCTION 

The resent development of ultrafast nonlinear opti- 
cal experiments 1 ' 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 * 12 ' 13 * 14 has stimulated 
immense interest in spintronics in semiconductors as it 
shows great potential of using the spin degree of freedom 
of electrons in place of/in addition to the charge degree of 
freedom for device application such as qubits, quantum 
memory devices, and spin transistors. 

In order to make use of the spin degree of 
freedom in semiconductor spintronics, it is cru- 
cial to have a thorough understanding of spin de- 
phasing mechanism. Three spin dephasing mecha- 
nisms have been proposed in semiconductors: 15 ' 16 the 
Ellit-Yafet (EY) mechanism, 17 ' 18 the D'yakonov-Perel' 
(DP) mechanism, 19 and the Bir-Aronov-Pikus (BAP) 
mechanism. 20 In the EY mechanism, the spin-orbit inter- 
action leads to mixing of wave functions of opposite spins. 
This mixing results in a nonzero electron spin flip due to 
impurity and phonon scattering. The DP mechanism is 
due to the spin-orbit interaction in crystals without in- 
version center, which results in spin state splitting of the 
conduction band at k =/= 0. This is equivalent to an effec- 
tive magnetic field acting on the spin, with its magnitude 
and orientation depending on k. Finally, the BAP mech- 
anism is originated from the mixing of heavy hole and 
light hole bands induces by spin-orbit coupling. Spin-flip 
(SF) scattering of electrons by holes due to the Coulomb 
interaction is therefore permitted, which gives rise to spin 
dephasing. The dephasing rates of these mechanisms fol- 
low polarized system are calculated in the framework of 
single particle approximation. 15 For GaAs, the EY mech- 



anism is less effective under most conditions, due to the 
large band gap and low scattering rate for high quality 
samples. The BAP mechanism is important for either 
p-doped or insulating GaAs. For n-doped samples, how- 
ever, as holes are rapidly recombined with electrons due 
to the presence of a large number of electrons, spin de- 
phasing due to the regular BAP mechanism is blocked. 
Therefore, the DP mechanism (or possibly the EY mech- 
anism under certain conditions) is the main mechanism 
of spin dephasing for n-type GaAs. 

It is important to note that all the mechanisms 
above either provide or are treated as SF scatter- 
ing. Spin-conserving (SC) scattering, such as the ordi- 
nary Coulomb scattering, electron-phonon and electron- 
nonmagnetic impurity scattering which has been exten- 
sively studied in connection with optical dephasing and 
relaxation, 21 are commonly believed to be unable to 
cause spin dephasing as the corresponding interaction 
Hamiltonians commute with the total spin operator. 

Recently Wu proposed a many-body kinetic theory 22 
to study the spin precession and spin dephasing in insu- 
lating ZnSe/Zni_ a; Cd a ;Se quantum well (QW), 23 n-typed 
bulk GaAs samples, 24 and n-typed GaAs (110) QW, 25 
under moderate magnetic fields in the Voigt configura- 
tion. Based on this many-body theory, he further showed 
that the SC scattering can also cause spin dephasing in 
the presence of inhomogeneous broadening. 22,24 ' 25 ' 26 ' 27 
This novel spin dephasing mechanism has long been over- 
looked in the literature. Differing from the earlier study 
of the spin dephasing which comes from SF scattering, 
the spin dephasing through inhomogeneous broadening 
is caused by irreversibly disrupting the phases between 
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spin dipoles and is therefore a many-body effect. 22,26 ' 27 
Very recently we have shown that this inhomogeneous 
broadening effect also plays an important role in the spin 
transport 28,29 . 

In this paper, we study the spin dephasing in n-doped 
GaAs (100) QW's. We calculate the SDT by numeri- 
cally solving the many-body kinetic equations with all 
the scattering included. Differing from the previous in- 
vestigation in the bulk case where we are only able to 
get the SDT qualitatively, 24 here we get the SDT quan- 
titatively thanks to the reduction of dimension in the 
momentum space. Moreover as we include the electron- 
electron Coulomb scattering in our calculation, for the 
first time we are able to study the spin dephasing with 
extra large initial spin polarization (up to 100 %) which 
has not been investigate both experimentally and theo- 
retically. 

We organize the paper as follows: We present our 
model and the kinetic equations in Sec. II. Then in Sec. 
III(A) we show the time evolution of the spin signal where 
we show the contribution of the Coulomb scattering to 
the spin dephasing. In Sec. III(B) we investigate how the 
SDT changes with the variation of the initial spin polar- 
ization. The temperature dependence of the SDT under 
different spin polarization is discussed in detail in Sec. 
III(C), where we also highlight the difference between 
the present many-body theory and the earlier simplified 
theory. In Sec. III(D) we show the magnetic field depen- 
dence of the SDT. Finally we discuss how the electron 
density affect the SDT. We present the conclusion and 
summary in Sec. IV. 



II. KINETIC EQUATIONS 

We start our investigation from an n-doped (1 00) GaAs 
QW with well width a. The growth direction is as- 
sumed to be z-axis. A moderate magnetic field B is 
applied along the x axis. Due to the confinement of 
the QW, the momentum states along z axis are quan- 
tized. Therefore the electron states are characterized by 
a subband index n and a two dimensional wave vector 
k = (k x , k y ) together with a spin index a. In the present 
paper, the subband separation is assumed to be large 
enough so that only the lowest subband is populated 
and the transition to the upper subbands is unimpor- 
tant. Therefore, one only needs to consider the lowest 
subband. For n-doped samples, spin dephasing mainly 
comes from the DP mechanism. 19 With the DP term in- 
cluded, the Hamiltonian of the electrons in the QW takes 
the form: 

H =J2{^+[9^B + h(k)}-^\clc^,+H I . (1) 



Here £k = k 2 /2m* is the energy of electron with wavevec- 
tor k and effective mass m*. a are the Pauli matrices. 
In QW system, the DP term is composed of the Dressel- 
haus term 30 and the Rashba term. 31,32 The Dresselhaus 
term is due to the lack of inversion symmetry in the zinc- 
blende crystal Brillouin zone and is sometimes referred to 
as bulk inversion asymmetry (BIA) term. For the (100) 
GaAs QW system, it can be written as 33 ' 34 



^ IA (k) 

C A (k) 
j, BIA 



= 7 ^((fc 2 )-fc 2 ), 
(k) = . 



(2) 



Here (fe 2 ) represents the average of the opera- 
tor — (J|) 2 over the electronic state of the low- 
est subband and is therefore (ir/a) 2 . 7 = 
(4/3)(m*/m cv )(l/^/2m*3E g )( V /^/l - ry/3) and rj = 
A/ (Eg + A), in which E g denotes the band gap; A rep- 
resents the spin-orbit splitting of the valence band; m* 
standing for the electron mass in GaAs; and Tfi^y is a 
constant close in magnitude to free electron mass mo- 16 
Whereas the Rashba term appears if the self-consistent 
potential within a QW is asymmetric along the growth 
direction and is therefore referred to as structure inver- 
sion asymmetry (SIA) contribution. Its scale is propor- 
tional to the interface electric field along the growth di- 
rection. For narrow band-gap semiconductors such as 
InAs, the Rashba term is the main spin-dephasing mech- 
anism; whereas in the wide band-gap semiconductors 
such as GaAs, the Dresselhaus term is dominant. In the 
present paper, we will take only the Dresselhaus term 
into consideration as we focus on the spin dephasing in 
GaAs QW. The interaction Hamiltonian Hi is composed 
of Coulomb interaction H ee , electron- phonon interaction 
H p h, as well as electron- impurity scattering Hi. Their 
expressions can be found in textbooks. 21 ' 35 

We construct the kinetic Bloch equations by the 
nonequilibrium Green function method 21 as follows: 



Pk,aa' — Pk,cra' l C oh Pk,crcr' Iscatt 



(3) 



Here /?k represents the single particle density matrix. 
The diagonal elements describe the electron distribu- 
tion functions pk,<T<r = /k<r- The off-diagonal elements 
Pk,|-i = Pk describe the inter-spin-band polarizations 
(coherence) of the spin coherence. 23 Note that p k _i 1 = 
p* i_i = p^. Therefore, /k±i and pk are the quantities 
to be determined from Bloch equations. 

The coherent part of the equation of motion for the 
electron distribution function is given by 
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-^"Icoh = - 2(j {[9^bB + h x (k)]Imp k + h y (k)Rep k } + 4erlm^ Vqp k+q p k , 



(4) 



where Vq = Aire 2 /[no(q + qo)] is the 2D Coulomb matrix element under static screening, qo = (e 2 m* / kq) ^ ct /k=o,cr 
and kq is the static dielectric constant. The first term on the right hand side (RHS) of Eq. (4) describes spin precession 
of electrons under the magnetic field B as well as the effective magnetic field h(k) due to the DP effect. The scattering 
terms of electron distribution functions in the Markovian limit are given by 



dfa,a 
dt 



SCatt = \ - 2?r 9QX S ( £ k - £k-q - ^Wa) [Nq qz x{fk a - /k-q CT ) + /k<x(l - /k-qa) - Rc(pkPk-q)] 

-2wNi Y U 2 S(e k - £ k -q) [jw(l - fa-qa) ~ Re(pkPk-q)] ~ 2lT ^ ^(ek-q ~ £k + £k' ~ £k' 
q qk'cr' 

(1 - /k-qo-)/k<r(l - /k'<7')/k'-qo-' + ^PkPk-qt/k'a' - /k'-q CT ') + Pk'Pk:'-q(/k-q<T - ,/kcr) 

-{k <-> k — q, k' <-> k' — q}, 



(5) 



r 



in which {k <-> k q, k' <-» k' — q} stands for the same 
terms as in the previous {} but with the interchange 
k <-» k q and k' <-> k' — q. The first term inside the 
braces on the RHS of Eq. (5) comes from the electron- 
phonon interaction. A stands for the different phonon 
modes, i.e., one longitude optical (LO) phonon mode, one 
longitudinal acoustic (AC) phonon mode due to the de- 
formation potential, and two AC modes due to the trans- 
verse piezoelectric field. gqq z x are the matrix elements of 
electron-phonon coupling for mode A. For LO phonons, 



=2 



{4™^/K/27^ 2 + q 2 z )]}\I{iq z )\ 2 with a = 

e2 \// i /(2^LO)( K °o 1 — Kq 1 )- Koc 1S the optical dielectric 
constant and ^ s the LO phonon frequency. The form 



factor \I(iq z )\ 2 = tt 2 sin 2 y / '[y 2 (y 2 - tt 2 ) 2 } with y = q z a/2. 
Aq 9i A = l/Icxp^q^AAsT 1 ) — 1] is the Bosc distribu- 
tion of phonon mode A at temperature T. The second 
term inside the braces on the RHS of Eq. (5) results 
from the electron-impurity scattering under the random 
phase approximation with Ni denoting the impurity con- 
centration. U 2 = {4^e 2 /[« ( g 2 + q 2 z )]} 2 \I {iq z )\ 2 
is the electron-impurity interaction matrix element with 
Zi stands for the charge number of the impurity. Zi is 
assumed to be 1 throughout our calculation. The third 
term is the contribution of the Coulomb interaction. Sim- 
ilarly, the coherent and the scattering parts of the spin 
coherence are given by 



<9p_k 

dt 



coh = I^PbB + ih x (k) + h y (k)](f w i -/ k _i)+i^Vq[(/ k+q i -/ k+q _i)pk-Pk+ q (/ki -/ k -i)L (6) 



dp ( 

-g^ Iscatt = \ n Yl Sqg^Ok - £k-q - ^qq z \) [Pk-q(/ k I + /k-±) + (/k-qi + /k-q-± - 2 )Pk - 27V q ^ A (/3k - Pk-q)] 
+nN i U*5(e k - £k- q ) [(/ k I + /k-|)Pk-q - (2 - / k _ q i - /k-q-l)Pk] 

q 

~H 7rV q 2(5 ( £k - < l _£k + £k ' " £ k'-q)( (/k-qiPk + Pk-q/ k -i)(/k'i -/k'-qi + /k'-i -/k'-q-i) 



qk' 

+Pk[(l - /k'i)/k-qi + (! - /k'-i)/k-q-i - 2Rc(p k ,p k -q)] - Pk-q[/ k 'i(l - /k'-qi) 

+ (1 - /k'-i)/k'- q -| - 2Rc(p k 'Pk'-q)] ) }-{k ~ k - q, k' ~ k' - q} . 



(7) 



The initial conditions are taken at t = as: /k<r|t=o = l/{ cx p[(£k — Pa)/kBT] + l} (9) 

Pk|t=o = (8) 
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where p, a is the chemical potential for spin a. The con- 
dition iii 7^ n_i gives rise to the imbalance of the clcc- 

2 2 

tron densities of the two spin bands. Eqs. (3) through 
(7) together with the initial conditions Eqs. (8) and (9) 
comprise the complete set of kinetic Bloch equations of 
our investigation. 



III. NUMERICAL RESULTS 



The kinetic Bloch equations form a set of nonlinear 
equations. All the unknowns to be solved appear in the 
scattering terms. Specifically, the electron distribution 
function is no longer a Fermi distribution because of the 
existence of the anisotropic DP term h(k). This term in 
the coherent part drives the electron distribution away 
from an isotropic Fermi distribution. The scattering term 
attempts to randomize electrons in k-space. Obviously, 
both the coherent part and the scattering terms have 
to be solved self-consistently to obtain the distribution 
function and the the spin coherence. 



2.05 





10.8 




12.9 


LOO 


35.4 meV 


* 

m 


0.067 m 


A 


0.341 eV 


E 9 


1.55 eV 


9 


0.44 







E 

% 1-95 



1.9 



1.85 



Y i i i i i i i i i i i i i i i 


i | i i i i - 
- 

- 


B 4 1 




Ni = 




T=120K 

, , , i , , , , i , , , , i 


, i , , , , 



10" 1 



-. IP 



-- 10- 



lO" 4 



20 40 60 
t(ps) 



80 100 



FIG. 1: Electron densities of up spin and down spin and the 
incoherently summed spin coherence p versus time t without 
taking account the Coulomb scattering for a GaAs QW with 
initial spin polarization P — 2.5%, at T = 120 K. Note the 
scale of the spin coherence is on the right side of the figure. 



A. Temporal evolution of the spin signal 



TABLE I: Parameters used in the numerical calculations 

We numerically solve the kinetic Bloch equations in 
such a self-consistent fashion to study the spin preces- 
sion between the spin-up and -down bands. We include 
electron-phonon scattering and the electron-electron in- 
teraction throughout our computation. As we concen- 
trate on the relatively high temperature regime (T > 
120 K) in the present study, for electron-phonon scatter- 
ing we only need to include clectron-LO phonon scat- 
tering. Electron-impurity scattering is sometimes ex- 
cluded. As discussed in the previous paper, 23 ' 36 irre- 
versible spin dephasing can be well defined by the slope of 
the envelope of the incoherently summed spin coherence 
p(t) = ^ k \pk\. The material parameters of GaAs for our 
calculation are tabulated in Table I. 37 The method of the 
numerical calculation has been laid out in detail in our 
previous paper on the DP mechanism in 3D systems. 24 
The difference is that here we are able to get the re- 
sults quantitatively in stead of only qualitatively as in 
our previous 3D case, thanks to the smaller dimension 
in the momentum space. Our main results are plotted 
in Figs. 1 to 10. In these calculations the total electron 
density N e and the applied magnetic field B are chose 
to be 4 x 10 11 cm~ 2 and 4 T respectively unless other- 
wise specified. The width of the quantum well is chosen 
typical to be 15 nm except in the last two figures. 



We first study the temporal evolution of the spin sig- 
nal in a GaAs QW at T = 200 K. In Fig. 1 we plot 
the electron densities in the spin-up and -down bands 
together with the incoherently summed spin coherence 
for Ni = 0. At t = 0, the initial spin polarization 
P = (N 1/2 - N_ 1/2 )/(N 1/2 + N_ 1/2 ) is 2.5%. In this 
calculation, the Coulomb scattering is not included. It is 
seen from the figure that excess electrons in the spin-up 
band start to flip to the spin-down band at t = due to 
the presence of the magnetic field and the DP term h(k). 
In the meantime the spin coherence p accumulates. At 
about 9.7 ps, the electron densities in the two spin bands 
become equal and the spin coherence reaches its maxi- 
mum. Then the spin coherence starts to feed back and 
the electron density in the spin-down band exceeds that 
in the spin-up band while p deceases. At about 18 ps, 
p reach its minimum, while the density difference in the 
two spin bands reaches its maximum again with the ex- 
cess electrons now in the spin-down band. Due to the 
the dephasing, the second peak is lower than the first 
one. This oscillation goes on until the amplitude of the 
oscillation becomes zero due to the dephasing. 

In Fig. 2, we plot the time evolution of electron den- 
sities in the two spin bands as well as the incoherently 
summed spin coherence for the same GaAs QW system 
as the previous one but taking the Coulomb scattering 
into account. The results without the Coulomb scatter- 
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FIG. 2: Electron densities of up spin and down spin and 
the incoherently summed spin coherence p versus time t with 
(solid curves) and without (dashed curves) taking account the 
Coulomb scattering for a GaAs QW initial spin polarization 
P — 2.5% at T = 120 K. Note the scale of the spin coherence 
is on the right side of the figure. 



ing are replotted as dashed curves in the figure. It is seen 
from the figure that for the first oscillation (t < 20 ps), 
the electron densities as well as the spin coherences are 
almost the same in the the presence and absence of the 
Coulomb scattering. As time goes on, the curves with 
the Coulomb scattering deviate from the ones without 
the Coulomb scattering. The decay rates of the excess 
spin density as well as the spin coherence are faster in 
the presence of Coulomb scattering. 

It is known that, the Coulomb scattering is important 
only when the electron distribution is divagated from the 
Fermi function. In the first few picoseconds, the electrons 
in the two spin bands flip to their opposite bands due to 
the magnetic fields. The buildup of the inhomogeneity of 
k in the electron distribution function comes from the DP 
term is marginal, and the electron distributions remain 
approximately the Fermi function. As time goes on, the 
effect of the DP term accumulates, the electron distri- 
bution functions divagate further and further away from 
the Fermi function and hence the Coulomb scattering be- 
comes more and more important. Consequently, the spin 
signal that with the Coulomb scattering differs from the 
one without the Coulomb scattering. It has been pointed 
out that the spin dephasing due to the DP term and 
the SC scattering comes from two effects. 24 The first one 
is widely discussed that is due to the anisotropic prop- 
erty of the DP term, which, combined with SC scattering 
gives rise to the effective SF scattering. 15 ' 19 In this case, 
inclusion of additional scattering enhances the momen- 



tum relaxation rate and consequently reduces the spin 
dephasing rate. 15,19 ' 22 It has been pointed out that the 
electron-electron Coulomb scattering, although does not 
contribute to the momentum relaxation rate, trends to 
reduce the spin dephasing 38 based on the similar analy- 
sis as in Refs. 15,19 and 22. The second is that the DP 
term itself also introduces an inhomogeneous broaden- 
ing, which in the presence of the SC scattering provides 
additional spin dephasing channel 22, 24,25,26,27 and there- 
fore results in a faster spin dephasing. Our calculation 
self-consistently solves the Kinetic Bloch equations and 
includes both effects. The result indicates that for the 
present condition the second effect is more important and 
therefore the combined effect by inclusion of the Coulomb 
scattering leads to the increase of the spin dephasing. 



B. Spin polarization dependence of the spin 
dephasing time 

We now turn to study the spin polarization dependence 
of the SDT. As our theory is a many-body theory and we 
include all the scattering, especially the Coulomb scatter- 
ing, in our calculation, we are able to calculate the SDT 
with large spin polarization. 

In Fig. 3, the SDT t is plotted against the initial spin 
polarization P for GaAs QW's with Ni = [Fig. 3(a)] 
and Ni = 0.1 N e [Fig. 3(b)] at different temperatures. 
The most striking feature of the impurity-free case is the 
huge increase of the SDT in low temperatures. For T = 
120 K, the SDT increases from 25 ps at low polarization 
to 720 ps at ~ 100 % polarization. In other words, the 
spin dephasing rate is decreased more than one order of 
magnitude when the spin polarization increases from 
to 100 %. It is also seen from the figure that the increase 
of the SDT is reduced with the increase of temperature. 
For T = 300 K, the SDT only gets an 80 % increase when 
the polarization increases from to 100 %. 

The dramatic increase in the t-P curve in the low tem- 
perature regime originates from the electron-electron in- 
teraction, specifically the Hartree-Fock (HF) self-energy 
[i.e., the last terms in the Eq. (4) and (6)]. As we know 
that the HF term itself does not contribute to the spin de- 
phasing directly. 22,27 However, it behaves as an effective 
magnetic field which can alert the motion of the electron 
spins and can therefore affect the spin dephasing by com- 
bining with the DP term. For small spin polarization as 
commonly discussed in the literature, the contribution of 
the HF term is marginal. However, when the polarization 
gets higher, the HF contribution becomes larger. For ex- 
ample, the magnitude of the effective magnetic field of 
HF term is larger than 40 T, which is about ten folds of 
the applied magnetic field, for the case when the tem- 
perature is 120 K and the spin polarization is ~ 100 %. 
Differing from the applied magnetic field which in the 
Voigt configuration only gets the transverse components 
(i.e. B x and B y ) and always causes the electrons to flip 
between spin-up and -down bands, the effective magnetic 
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FIG. 3: Spin dephasing time r versus the initial spin po- 
larization with different impurity concentration and different 
temperatures. The impurity densities in (a) and (b) are and 
0.1 N e respectively. The lines are plotted for the aid of eyes. 



field formed by the HF term contains a longitudinal com- 
ponent [Bp-^(k)] which can effectively reduce the "de- 
tuning" of the spin-up and -down electrons, and thus 
strongly reduces the spin dephasing. In order to show 
this detuning effect, we remove the longitudinal compo- 
nent of the effective magnetic field Bp (k) and recalcu- 
late the SDT at the temperature of 120 K for different 
initial spin polarization for an impurity free sample. The 
result is plotted in Fig. 4. From the figure one can see 



FIG. 4: Spin dephasing time r versus the initial spin polar- 
ization P for at T = 120 K and Ni = 0. Circle (•): With the 
longitudinal component of HF term included; Diamond (♦): 
Without the longitudinal component of HF term include. The 
lines are plotted for the aid of eyes. 

that when the longitudinal component of the HF term 
is removed, the dramatic increase in the t-P curve dis- 
appears and the SDT is insensitive to the initial spin 
polarization. 

When the temperature increases, for a given initial spin 
polarization the HF term becomes smaller as the elec- 
trons are distributed to a wider range in the k-space. 
Therefore the effect of the HF term becomes smaller 
too. Consequently the increase in the t-P curve becomes 
slower. 

The t-P curve gets flatter when the impurities are in- 
troduced. It is shown from Fig. 3(b) that, when the 
density of impurity is large, say Ni =0.1N e , the fast rise 
in the t-P curve in low temperature regime still remains. 
Nevertheless the rate of increase is much smaller than the 
corresponding one when the impurities are absent. 

To further reveal the contribution of the impurity to 
the dephasing under different conditions, we plot the 
SDT as a function of the polarization for different im- 
purity levels at T = 120 K in Fig. 5. The figure clearly 
shows that when the impurity concentration increases, 
the slope of the t-P curve becomes smaller. This is be- 
cause that impurity scattering reduces the HF term and 
the effect of the longitudinal component of the HF term 
is also reduced. Consequently, the increase of SDT with 
the polarization is reduced. 

It is interesting to note that contrary to the high polar- 
ization regime where the SDT decreases with the impu- 
rity concentration, the spin dephasing is reduced by the 
impurity scattering in the low polarization regime. As 
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FIG. 5: Spin dephasing time r versus the initial spin polariza- 
tion P at T — 120 K for different impurity levels. Circle (•): 
Ni = 0; Diamond (♦): iV 4 = 0.01iV e ; Square (■): N t = 0.1iV e . 
The lines are plotted for the aid of eyes. 



we pointed out before that the impurities affect the spin 
dephasing in two ways. 24 On the one hand, the electron- 
impurity scattering provides a new spin dephasing chan- 
nel by combining with the DP term 19 ' 24 to give an effec- 
tive SF scattering through the inhomogeneous broaden- 
ing introduced by the DP term. This gives rise to the 
enhancement of the dephasing. On the other hand, the 
scattering also redistributes the electrons in the momen- 
tum space and leads them to an isotropic distribution. 
Therefore, the scattering can suppress the anisotropy 
caused by the DP term, consequently the effective SF 
scattering. Moreover, the suppression of the anisotropy 
also corresponds to the reduction of the inhomogeneous 
broadening. Both lead to a smaller spin dephasing. Our 
results indicate that in the low temperature and the low 
polarization regime, the impurities tend to reduce the 
spin dephasing. 



C. The temperature dependence of the spin 
dephasing time 

Above we discussed the dependence of spin dephas- 
ing on initial spin polarization for different temperatures. 
Now we turn to the temperature dependence of the SDT 
under different initial spin polarizations. From Fig. 3(a) 
and (b) in Sec. III(B), one can see that for small spin 
polarization, the SDT increases with the temperature. 
Whereas in high polarized regime, the SDT decreases 
with the temperature. For moderate polarization, the 



temperature dependence is too complicated to be de- 
scribed by a monotonic function of temperature. Un- 
der certain condition, the SDT can be insensitive to the 
temperature, e.g. the SDT is almost unchanged with 
the temperature when the polarization P = 75 % and 
N { = 0.1N e . 

To see more detail of how the spin dephasing depends 
on the temperature, we replot in Fig. 6(a) and (b) the 
SDT shown in Figs. 3 and 4 as a function of the tem- 
perature for different impurity levels and different spin 
polarizations. It is seen from the figure that, for low 
polarization, the SDT increases systematically with the 
temperature for all impurity levels. This property is op- 
posite to the results of earlier simplified treatments of 
the DP effect, where it was predicted that the spin life- 
time decreases with the increase of temperature in the 
2D system. 39 - 40 The SDT based on the simplified model 
is given by 26,39 ' 41 



1 = f~dE k (f k i-f k _i)T(k) 
f °°dE k (f kl -f k _ h ) 



(10) 



in which 



T(k) 



2r 1 (fc)[( 7 (fc 2 )) 2 fc 2 -i 7 (fc 2 )fc 4 



+ 



16 



(11) 



and 



a{E k ,6)[l - cos(n6)]d6 . (12) 



a(E k ,6) stands for scattering cross-section. For compari- 
son, we plot the SDT predicated by the earlier model and 
by our present many-body theory in the inset of Fig. 6(a). 
From the inset one can see that the SDT predicated by 
the earlier model is about one order of magnitude larger 
than the one predicated by our theory. In the mean time, 
the SDT of the earlier mode drops dramatically with the 
increase of the temperature. Nevertheless, in our many- 
body treatment, it rises slightly with the temperature. 

The recent experiments show that the SDT in n- 
type quantum wells increases slightly with the increase 
of temperature, 42 or is almost unchanged with the 
temperature. 43 It is also noted that the SDT's measured 
in the experiment 43 for GaAs QW are with the regime of 
values predicted by our theory, but one order of magni- 
tude smaller than those by the earlier model. Moreover, 
it has been reported very recently that there is a big dis- 
crepancy between the earlier simplified treatment of the 
DP effect and the experiment on the SDT. It is shown 
that the theoretical SDT is about one order of magnitude 
larger than the experiment data in n-type bulk GaAs for 
certain electron densities. 44 

It is seen that our many-body result is better than 
the earlier simplified treatments both qualitatively and 
quantitatively. The reason that our model is more pre- 
cise than the earlier one lies on the fact that the earlier 
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model is based on the single particle picture which does 
not count for the dephasing due to the inhomogeneous 
broadening inherited in the DP term, which is the re- 
sult of the many body effect. 22 ' 24,25 ' 26,27 By comparing 
the theoretical SDT predicated by the two models, we 
can see that the spin dephasing due to the inhomoge- 
neous broadening is much more important. In the case 
we calculated, the spin dephasing is dominated by the 
inhomogeneous broadening. Therefore, it is easy to un- 
derstand why the earlier simplified treatment of the DP 
mechanism gives much slower spin dephasing. 

The temperature dependence of the SDT can be under- 
stood once the spin dephasing due to the inhomogeneous 
broadening is taken into account: When the tempera- 
ture increases, the inhomogeneous broadening is reduced 
as the electrons are distributed to the wider fc-states. As 
a result, the number of electron occupation on each k 
state is reduced. It is further noted that this reduction 
is mild as a function of the temperature. Therefore, the 
temperature dependence is quit mild unless it is within 
the regime where the HF term plays an important role 
in the spin dephasing. 

In the region where the HF term is important, in addi- 
tion to the above mentioned two effects of the tempera- 
ture on the spin dephasing, the temperature dependence 
of the HF term should also be taken into account. In 
general, the temperature dependence of the SDT due to 
the combination of these three effects is too complicated 
to be described by a monotonic function. We replot the 
SDT as a function of the temperature in Fig. 6(b), for a 
typical high polarization P = 75 %. We can see that due 
to the reduction of the HF term, the reduce of the de- 
tuning is removed and SDT drops dramatically with the 
increase of the temperature in the impurities free sam- 
ple. While for the system with impurity concentration 
Ni = 0.1N e , the HF term is not as important as in the 
impurities free sample, and the SDT is insensitive to the 
temperature. 



D. Magnetic field dependence of the spin 
dephasing 
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We now investigate the magnetic field dependence of 
the spin dephasing. In Fig. 7, we plot the SDT versus 
the applied magnetic field for different impurity levels 
and different spin polarizations. It is seen from the fig- 
ure that in small polarization regime for both impurity 
free and doped samples, the SDT is almost a constant 
when the magnetic field varies from 1 T to 8 T. Whereas, 
in the high polarization regime, the SDT increases with 
the magnetic field when the magnetic field varies from 
1 T to 4 T and then saturates when the magnetic field 
is larger than 4 T. It is noted that the transverse mag- 
netic field imposes two effects on the spin dephasing. One 
is that in the presence of a magnetic field, the electron 
spins undergo a Larmor precession around the magnetic 
field. This precession suppresses the precession about the 



FIG. 6: Spin dephasing time r versus the temperature T with 
spin polarization P — 2.5 % (a) and P — 75 % (b) under 
two different impurity levels. Circle (•): Ni = 0; Square 
(■): Ni = 0.1iV e . The lines are plotted for eye aid. The 
SDT predicated by the simplified treatment of DP term (solid 
curve) and our model (circle) for Ni = is plotted in the inset 
of (a) for comparison. 



effective magnetic field h(k) of the DP term. 15 ' 45 There- 
fore the SDT increases with the magnetic field. How- 
ever, in the presence of the transverse magnetic field, 
the spin precession frequency between the spin-up and 

-down band is \J {c^ibB — h x (h)) 2 + /i 2 (k). Hence the 
transverse magnetic field may introduce an additional 
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FIG. 7: Spin dcphasing time r versus the applied magnetic 
field for different spin polarizations and different impurity 
levels. Solid curve with dots: N = 0, P = 2.5 %; Solid 
curve with squares: Ni = 0.1N e , P = 2.5 %; Dashed curve 
with dots: iV» = 0, P = 75 %; Dashed curve with squares: 
Ni = O.liVe, P = 75 %. 

inhomogeneous broadening for the electrons and conse- 
quently results in a shorter SDT. In the cases we study, 
the SDT is almost unchanged with the magnetic field 
with the combined contributions from these two effects. 

In additional to the above mentioned effects of the 
magnetic field on the spin dephasing, one can further see 
from Fig. 7, that for large spin polarization, the magnetic 
field enhances the SDT. As we mentioned before, for large 
spin polarization, the contribution from the HF term is 
important. Therefore, the magnetic field can further af- 
fect the spin dephasing rate through the enhancement of 
the HF term in the high initial spin polarization regime. 
It is shown in the figure that, for high initial spin polar- 
ization, when the magnetic field increases, the enhanced 
HF term results in a fast decrease of spin dephasing rate. 
Nevertheless, after the magnetic field is high enough, the 
HF term arrives it maximum. As a result, the SDT sat- 
urates with the magnetic field. 

E. Electron density dependence of spin dephasing 

We now turn to the electron density dependence of the 
spin dephasing. In Fig. 8 we plot the SDT as a func- 
tion of the electron density for both low and high spin 
polarizations. It can be seen from the figure that for 
low polarization, the SDT decreases with electron den- 
sity: from 70 ps for N e = 4 x I0 11 cm~ 2 to 4 ps for 
N e = 40 x 10 11 cm~ 2 . This is because with the increase 



FIG. 8: Spin dephasing time r versus the total electron den- 
sity N e for a GaAs QW with T = 200 K, B = 4 T, N = 
and P = 2.5 % (•) and P = 75 % (♦). The lines are plotted 
for the aid of eyes. 



of the electron density, more electrons are distributed 
at large momentum states and strengthen the DP ef- 
fect as the DP term increases with the momentum. For 
the high spin polarization, the SDT also decreases with 
the electron density, however, with a much faster speed. 
When the initial spin polarization is 75%, the SDT de- 
creases from 280 ps for N e = 4 x 10 11 cm~ 2 to 8 ps for 
A^e = 40 x 10 11 cm~ 2 . This is understood that when the 
electron density increases, the effect of HF term becomes 
less important comparing to the increase of the DP term. 



F. The effect of quantum well width on the spin 
dephasing 

It is noted from the experiments that the SDT of quan- 
tum wells is much smaller than that of bulk material. 
This is because the DP term in the quasi-two dimensional 
quantum well contains a term which is proportional to 
the (fc 2 ) = (n /a) 2 . For a regular quantum well in the 
order of 10 nm, this term is much larger than the square 
of Fermi vector hp, when the electron density is up to the 
order of 10 11 cm~ 2 . Therefore the DP term in the quan- 
tum well is greatly enhanced. Moreover, the smaller the 
width is, the larger the DP term becomes and the faster 
the SDT turns to be. In this subsection we study the 
effect of the quantum well width on the spin dephasing. 
In Fig. 9 the SDT is plotted as a function of the initial 
spin polarization for three different widths. Moreover, 
in Figs. 10(a) and (b) we plot the SDT as a function 
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of the temperature for two initial spin polarizations in 
three different quantum wells. These figures clearly show 
that, for all situations, the spin dephasing rate decreases 
with the increase of the well width. Especially when the 
width is decreased by only 50 % from 20 nm, the SDT is 
reduced by one order of magnitude. 



1500 



1200 




FIG. 9: Spin dephasing time r versus the initial spin polar- 
ization P at T = 120 K and Ni — in three quantum wells 
with different widths: Diamond (♦): a — 10 nm; Circle (•): 
a — 15 nm; Square (■): a — 20 nm. The lines are plotted for 
the aid of eyes. 



IV. CONCLUSION 

In conclusion, we have performed a systematic investi- 
gation of the DP effect on the spin dephasing of n-typed 
GaAs QW's for high temperatures under moderate mag- 
netic fields in Voigt configuration. Based on the noncqui- 
librium Green's function theory, we derived a set of ki- 
netic Bloch equations for a two-spin-band model. This 
model includes the electron-phonon, electron-impurity 
scattering as well as the electron-electron interaction. By 
numerically solving the kinetic Bloch equations, we study 
the time evolution of electron densities in each spin band 
and the spin coherence - the correlation between spin-up 
and -down bands. The spin dephasing time is calculated 
from the slope of the envelope of the time evolution of 
the incoherently summed spin coherence. We therefore 
are able to study in detail how this dephasing time is af- 
fected by spin polarization, temperature, impurity level, 
magnetic field and electron density. Different from the 
earlier studies on spin dephasing based on the single par- 
ticle model which only considers the effective SF scatter- 
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FIG. 10: Spin dephasing time r versus temperature T in 
three quantum wells with different widths: Diamond (♦): 
a = 10 nm; Circle (•): a — 15 nm; Square (■): a = 20 nm. 

(a), the initial spin polarization P — 2.5 %; (b) P — 75 %. 
The lines are plotted for the aid of eyes. 



ing, our theory also takes account of the contribution of 
many-body effect on the spin dephasing. 22,24,25,26,27 j n 
fact, for the n-typed semiconductors and the spin polar- 
ization studied in the experiments, this many-body de- 
phasing effect is even more important than the effective 
SF scattering as it can be one order of magnitude larger 
than the later. Equally remarkable is that, as we include 
all the scattering, especially the Coulomb scattering in 
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our many-body theory, now we are able to calculate the 
spin dcphasing with extra large (up to 100 %) initial spin 
polarization. 

It is discovered that the SDT increases dramatically 
with the initial spin polarization. We stress here that 
the SDT is defined to be the inverse of the spin dephas- 
ing rate, instead of the total life time which naturally 
increases with the initial spin polarization. At low im- 
purity level and for low temperature, the magnitude of 
SDT increases more than one order when the initial po- 
larization goes from to 100 %. It is discovered that this 
dramatic increase originates from the HF contribution 
of the electron-electron Coulomb interaction. The HF 
term forms an effective magnetic field which affects the 
spin dephasing. As the longitudinal component of the 
HF term effectively reduces the "detuning" between the 
spin-up and -down bands, the spin dephasing becomes 
much slower in high polarization region. Due to the fact 
that the dramatic increase comes from the HF term, the 
magnitude of the incrcasemcnt is therefore affected by all 
the factors that influences the HF term, such as temper- 
ature, impurity scattering, magnetic field as well as the 
electron density. When the temperature or the impu- 
rity density increases, the effective magnetic field formed 
by the HF term for a given initial spin polarization de- 
creases. Therefore, the increase in t-P curve is reduced, 
and the t-P curve becomes flatter. Nevertheless when 
the magnetic field increases, the HF term is enhanced 
and the increasement in the t-P is enhanced accordingly. 
When the magnetic field increases beyond 4 T, the HF 
saturates. Consequently after 4 T, the increase in t-P 
curves only gets slight changes with the magnetic field. 

For low spin polarized regime, the SDT increases with 
the temperature. This is contrary to the result of ear- 
lier simplified single-particle calculation where the SDT 
always decreases with the increase of the temperature. 
Moreover, the SDT predicted by our many-body calcula- 
tion is one order of magnitude faster than the earlier re- 
sult. We show that our results are in agreement with the 
experiments both qualitatively and quantitatively. The 
physics of this feature is due to the additional many-body 
spin dephasing channel due to the inhomogeneous broad- 
ening provided by the DP term, which by combining with 
the SC scattering also causes spin dephasing. In the sit- 
uation we studied, the spin dephasing is dominated by 
the many-body dephasing effect. With the increase of 
the temperature, the inhomogeneous broadening reduces 
and the SDT increases. 

The effect of the electron-impurity scattering on the 
spin dephasing is also studied. The SDT increases in 
the presence of the impurity scattering. This is because 



the impurity scattering redistributes the electrons to an 
isotropic state which in turn reduces the spin dephasing. 

In high spin polarization region, in additional to the 
above mentioned role the scattering plays in the spin de- 
phasing, the scattering also affects the HF term. This 
brings more complication in the study of the spin de- 
phasing. In general the SDT can not be described by a 
monotonic function of the temperature or the impurity 
concentration in high polarization region. For a typically 
high initial spin polarization, P = 75 %, the SDT de- 
creases dramatically with the temperature as the reduc- 
tion of the detuning is suppressed when the temperature 
increases. Whereas when the impurity concentration is 
0.1A^ e , the SDT is insensitive to the temperature. 

As the magnetic field causes the electron spins to pre- 
cess about it, this precession will suppress the precession 
about the effective magnetic field h(k) originated from 
the DP effect. However, the transverse magnetic field in 
the Voigt configuration also introduces additional inho- 
mogeneous broadening in the momentum space and en- 
hances the spin dephasing. As a result of the combining 
effects, for the condition we studied, the spin dephasing 
is almost unchanged with the magnetic field in low spin 
polarization region. The magnetic field also enhances the 
HF term, therefore the t-B curve gets a faster increase 
in the high polarization region when the magnetic field 
is smaller than 4 T. In the region of the magnetic field 
higher than 4 T, the HF term achieves its maximum. The 
effect of the magnetic field on the HF term saturates, the 
SDT saturates accordingly. 

Our calculation also shows that when then electron 
density increases, the SDT decreases. This is because 
with the increase of the electron density, more electrons 
are distributed at larger momentum states and conse- 
quently the DP term is strengthened. While as the width 
of the quantum well increases, the SDT decreases as the 
DP term reduces. 

In summary we have performed a thorough investiga- 
tion of the spin dephasing in n-typed GaAs QW's for 
high temperatures. Many new features which have not 
been investigated both theoretically and experimentally 
are predicted in a wide range of parameters. 
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